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' The heavy baryon Aq (Q = 6 or c) can be regarded as composed of a heavy quark and a scalar hght 

diquark which has good spin and isospin quantum numbers. In this picture we estabhsh the Bethe- 
^ Salpeter (BS) equation for Aq to second order in the l/mg expansion. With the kernel containing 

both the scalar confinement and the one-gluon-exchange terms we solve the BS equation numerically. 
' The value of the spin-dependant form factor for the matrix element {Ab{v, s)\b'y'^'y5b\Ab{v, s)), €b, 

which is non-zero at order is obtained numerically from our model. 
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I. Introduction 

In the past decades there has been much progress in heavy flavor physics due to the discovery of the new 



> 
00 

flavor and spin symmetries SU{2)f x SU{2)s in the heavy quark hmit and the establishment of the heavy 
quark effective theory (HQET) The Large Hadron CoUider (LHC) will provide much more data for heavy 



hadrons, and hence it will be able to test the standard model (SM) more accurately. One may expect more 



> precise measurement of Cabibbo-Kobayashi-Maskawa matrix elements such as Vub in the near future. 



There have been extensive studies in literature on inclusive semileptonic decays of bottom hadrons, iJf, — > 



Xeve {X represents all hadrons in the final states) [2|-[l3|, especially since the establishment of HQET. These 



studies include corrections to the leading order results both from perturbative QCD (as(mf,)) terms and from 
nonperturbative terms which are suppressed by powers of mf,. It was pointed out that there are no 1/mb 
corrections to leading order in the 1 / mf, expansion for the differential decay widths of inclusive semileptonic 
decays of bottom hadrons, dT / dq^dE^, where q is the total momentum of the electron and the neutrino and ^-eis 
the electron energy in the A;, 's rest frame Q. Bigietal. studied corrections to the decay width dF/dE'e Q]- 

Manohar and Wise analyzed extensively corrections to dV / dq^dE^ for unpolarized bottom hadron i/f, and 

for polarized A^ In recent years, theoretical calculations for the inclusive semileptonic decay widths and for 
the moments of inclusive observables have been carried out to order 1/m^ and a^/?o (/?o — 11 — 2n//3, nj is 
the number of quark flavors) 

At leading order in as{mb) the polarized differential semileptonic decay rate for Af, XucDe (where Xu 
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represents all hadrons containing an up quark) is 



dq^ dEe dEy c? cos 



27r3 



Wiq^ + W2{2E,E, - + W:iq^{E, - E,) 



47^3 



-G6(g" - 4SeK) - Ggg" - Gm { Ee - E, 



Giq^ + G2 ( 2EeE, - ) + Gag' | (£;e - E,)\ 



2E, 



Ep, + 



(1) 



where the electron mass is neglected, j = u,c, is the angle between the electron three-momentum and the 
spin vector of A;, in the A^'s rest frame, E^ is the neutrino energy in the rest frame of A;,, q^ = (pe + p^Y is 
the invariant mass of the lepton pair, the kinematic variables are to be integrated over the region q^ < AE^E^, 
Wi, W2, and W3 are the form factors in spin-independent hadronic tensor W'^'^, Gi, G2, G3, Gg, Gg, and Gg 
are the form factors in spin-dependent hadronic tensor Wli'^ . W^"^ and Wg'^ have the following forms 



and 



= -q-s[-g^'^Gi+v^'v''G2~ie^''"^v^qpG^ + q^q''Gi + {q^v'' + q''v^)G^] 



(2) 



(3) 



where is the spin vector of Ah and is the velocity of A;,. 

When the electron mass is ignored, the form factors W4, W5, G4, G5, and G7 in Eqs.(2)(3) do not contribute to 
the differential inclusive semleptomic decay rate of polarized Af, [8]. The expressions for other form factors, Wi, 
W2 , W3 , Gi , G2 , G3 , Ge , Gs, and Gg , were obtained to second order in the 1 /m?, expansion in Ref. Q . It was 
shown that l/m^ corrections were characterized by two parameters, fi^ and £{,, for polarized Af, decays (the 
spin energy is zero for Ah), fi^ is the kinetic energy which is defined as 



{Ab\K{iD^fhy\Ab) 
2mii 



(4) 



where denotes the field of b quark in HQET, D'^ = I?^ — v^^v ■ D with being the covariant derivative. 
The value of was calculated for the first time in Ref. [14| . e;, is defined as 



{Ai,{v,s)\bj^'-f5b\Ab{v,s)) ^ {1 + eb)uA,{v,s)j^'-f5UA^{v,s) 



(5) 



where u\^{v,s) is the Dirac spinor of Af, with helicity s. ef, is equal to zero at leading order in the l/mf, 
expansion due to the heavy quark symmetry. It is also zero at the first order in 1 / mi, expansion because of the 



current conservation [15|. It is the aim of the present paper to calculate the value of et, at order l/m^. This 
will be important for the precise measurement of Vub in polarized Af, decays. 

When the quark mass is very heavy compared with the QCD scale Aqcd, the light degrees of freedom in 
a heavy baryon Aq {Q = b or c) become blind to the flavor and spin quantum numbers of the heavy quark 
because of the SU{2)f x SU{2)s symmetries. Therefore, the angular momentum and flavor quantum numbers 
of the light degrees of freedom become good quantum numbers which can be used to classify heavy baryons, 
and Aq corresponds to the state in which the angular momentum of the light degrees of freedom is zero. So 
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it is natural to regard the heavy baryon as composed of one heavy quark and a hght diqiiark. When l/mg 
corrections arc taken into account, since the isospin of Aq is zero, the isospin of the light degrees of freedom is 
also zero. Therefore, the spin of the light degrees of freedom should also be zero in order to guarantee that the 
total wave function of Aq is anti-symmetric. Hence the spin and isospin of the light degrees of freedom are still 
fixed even when l/mg corrections are taken in account. Therefore, we still treat Aq as composed of a heavy 
quark and a light diquark. 

Based on the above picture, the three-body system is simplified to a two-body system. We will establish 
the Bcthe-Salpetcr (BS) equation for Aq in this picture to order 1/m.Q. Then, we will solve this equation 
numerically by assuming that the kernel contains the scalar confinement and the one-gluon-exchange terms. 
Furthermore, we will express in terms of the BS amplitude of A^ in order to calculate the value of to 
second order in the l/nib expansion. 

The remainder of this paper is organized as follows. In Sec. II we will establish the BS equation for Aq 
to second order in the 1/toq expansion and discuss the form of its kernel. Then, we will solve the equation 
numerically. In Sec. Ill we will express £(, in terms of the BS amplitude of A;, and calculate the value of €b for 
different model parameters. Finally, Sec. IV is reserved for summary and discussion. 



II. Bethe-Salpeter equation for Aq to second order in the l/mg expansion 

As we discussed in Introduction, Aq is regarded as a bound state of a heavy quark and a light scalar diquark. 
Hence we can define the BS amplitude of Aq as the following: 

X{xux2, P) = {0\Tij{xr)<pix2)\AQ{P)), (6) 

where ip{xi) and ip{x2) are the field operators of the heavy quark at position xi and the light scalar diquark at 
position X2, respectively, P = uiaqV is the momentum of Aq, and v is its velocity. Let rriQ and mo represent the 
masses of the heavy quark and the light diquark in the baryon, Ai = , X2 = m^+mo ' ^^'^ ^ represents 

the relative momentum of the two constituents. Then, the BS amplitude in momentum space is defined as 

x(xi,X2,P) =e^^^|^e^^-XpW, (7) 

where X = Xixi + X2X2 is the coordinate of the center of mass and x = xi — X2- 

It is straightforward to prove that the BS equation for Aq has the following form in momentum space: 

Xpip) = SriPi) J ^K{P,p,q)xp{q)SD{p2), (8) 

where pi = XiP + p, p2 = —X2P +p are the momenta of the heavy quark and the light scalar diquark, 
respectively, K{P,p,q) is the kernel which is defined as the sum of two particle irreducible diagrams, Sf{pi) 
and Sd{P2) are propagators of the heavy quark with momentum p\ and the light diquark with momentum p2, 

Sf{pi) = (9) 

- niQ + le 

Sd{P2) = -2 V— ^. (10) 

P2- mjj + ze 
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In order to solve the BS equation for Aq to second order in the l/mg expansion, we rewrite Eq. (8) as the 
following: 



Xop(p) H Xip(p) H tX2p{p) = [Sof{pi)^ Sif{pi) A 2-'S'2f(pi) / ■j^—-j[KQ{P,p,q) 

mq m^Q \ tuq ui^q J J (27r)4 ^ 

+ — KiiP,p,q) + ^K2iP,p,q)) X (xop{q) + —Xip{q) 

+ ^X2Piq))SD{p2), (11) 

where we have expanded the BS amplitude, the propagator of the heavy quark, and the kernel to l/mg. 

It is easy to show that Sd remains unchanged in the 1/toq expansion. Then, by comparing the two sides of 
Eq. (11) at each order in l/mg, we have following equations: 
to leading order in the 1/toq expansion: 

j^--^Ko{P,p,q)xop{q)SDiP2), (12) 



to first order in the l/mq expansion: 



/d 
j^^--^KQ{P,p,q)xop{q)SDiP2) + Sqf{pi) j -p^Ki{P,p,q)xQpiq)SDiP2) 

d^q 



d^q 
(2^ 



+5*0^(^1) J j^;^Ko{P,p,q)xipiq)SD{p2), (13) 

to second order in the 1 / mq expansion: 

/d^ f d^ 

-j^-^Ko{P,p,q)xop{q)SD{p2) + Sif{pi) j -j^^-^Ki{P,p,q)xop{q)SD{p2) 

j^^-^Ka{P,p,q)xip{q)SD{P2) + Sqf{pi) j j^^K2{P,p,q)xopiq)SD{P2) 



d^ f d^ 

-SoFiPi) J j^^^Ki{P,p,q)xip{q)SD{P2) + Sqf{pi) j j^^Kq{P,p, q)x2piq)SD{P2)- (14) 

In the following we will use the variables pi = v ■ p — A2TOaq, Pt — P ~ {v ■ p)v. Then, in leading order, first 
order and second order in the 1 /mg expansion, we have the following equations for the heavy quark propagator 
using the relation pi = AiP +p 



mq expa 

m m 



sopip.)- ^, j':^^ (15) 

2[pi + Eo + mo + le) 

9 rr,^ = -fci)+^ »(2gl-P^)(l+?<) 

4 ^ 2{pi + Eo + mD + ie) A{pi + Eo + mo + ie)^' ^ ' 



Q I \ i{pi+Eo + mD){l + i) ift iPt 

J2F[Pi) = 



8 4 A{pi + Eq + mo + ie) 

2iE2{l + /') + i{2Ei - pl)f^ i(2Ei ~ piY{l + f) 



A{pi + So + + ieY 8{pi + Eo + mD+ ie)^ ' 
where Eq, Ei, and E2 represent the binding energies, which satisfy the following relation: 



(17) 



mAo = wq + + £^0 + — — + — (18) 
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Sd{P2) remains unchanged in the l/mg expansion, and we rewrite it as the following using the relation 

P2 = -X2P + P- 



Sd{p2) = 



pf -W^ +ie' 



(19) 



where Wp — vTt 





[1, a 

FIG. 1: Diquark-gluon-diquark vetex, and a are the Lorentz and color indices of the gluon, respectively. 



We assume the kernel has the following form: 

-iK^^l® IVi + ® (j)2+P2YV2, 
-iKi = l® 1^3 + 7^ ® {P2+p'2YVi, 



(20) 
(21) 



iK2 



1 ® IV5 +7'' ® {,P2+p'2YVq, 



(22) 



where the first terms on the right hand sides of Eqs. (20), (21), and (22) arise from the scalar confinement 
and the second ones are from the one-gluon-exchange diagram, p2 and p'2 are the momenta of the light diquark 
attached to the gluon. and 7'' in Eqs. (20)-(22) are from the vertex of the heavy quark and the gluon. It 
is noted that 7^ in Eg. (20) becomes because of heavy quark symmetry (see FIG. 1). The forms of Vi and 
V2 were given before [la|- In this work, we assume that V3 and 14, V5 and Ve have the same forms as Vi and 
V2, respectively. However, V3 and V4, V5 and Vis are suppressed with respect to Vi and V2 by order Aqcd/itiq 
and Agp^/mQ, respectively. Note that in Eq. (20) is replaced by 7^ in Eqs. (21) and (22) since there is no 
heavy quark spin symmetry when we include l/mg corrections. 

Combining Eqs. (12)-(22), we have the BS equations for Aq in leading order, first order, and second order in 
the l/mg expansion: 



Xop(p) = 



2{pi + niD + E(i + is) J (27r) 



i[i (g> iVi +Vf,(s {P2+P2TV2]xop{q) 



pi-Wl+ie 



(23) 
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xip(p) = 



2{pi + mo + Eo + le) J (27r)4 



i[l ® IV3 + 7m ® (P2 + P2)^^4]xop(g) 



- W^l + ie 



2{pi + mD+Eo + ie) 

4:{pt + mo + Eq + ief 
dS 



(27r)4 
1 - 



i[i ® iVi +Vf,® {p2 + P2TV2]xip{q) 



2{pi +mD + Ea+ ie) 



X I 



(2^)^ 



1(g) IVi + v^, (g) ip2+P2)'"V2]xop{q) 



pf-W^ + ie' 



(24) 



X2P{P) = 



l{l+f) 



dS 



i[l ® IVl+V^, (E) {P2 + P2)'"V2]X2p{q) 2_„.2 , . 



2(p( + rriD + iio + «e) J (27r)^ pf - Wp + 



d'^q 



+1 



2{pi +mD+Eo+ le) J (27r)4 

(-2i;i+p^)(i + ^) 1-^ 



i[l ® 1^5 + 7m (P2 +P2)^^6]X0p(g) 



2(Pi +mD+Eo+ ie) 



i[l eg) IFi +Vf,(E) ip2 + P2)'"V2]xipiq) 



+i 



(-2i?i+p2)(i + ^) 1 



pf -W^ + ie 



A{pi +mD+EQ + ie)2 4 



2(p; + mz3 + £0 + «e) 



i[l <E) IV3 +7m ® (P2 +P2)^^4]xop('7)- 



(2^)4"L^-"- ' '^^^"^ ' '"^"""^"^pf-W^I 



+i 



{pi + niD + Eo){l + ■ 



{2E^-pf)f^ 



-P} 



2E2{l + i!) 



4 4(p; + mo + Eo + ie) A{pi + mo + Ea + ie)^ 
{2E,~pf)^{l + f) ] f dS 



Mpi + mo + Eo + ie)^ 8{pi + mo + Eq + ie)^ 
X Xop(g) 



(27r; 



^i[l(E) IVi + Vf, (g) (P2 + P2)^V2] 



Pi - ^p + 



(25) 



It is noted that Eqs. (23) and (24) arc just the equations obtained before in the heavy quark limit [16;], and 
to first order in the 1/toq expansion 17]. 
In general, xp{p) can be expanded as 



Xp{p) = {A + Bi, + Cf + Dff)uAQ {v, s), 
where A, B, C, D are Lorentz scalar functions. It is easy to prove that 

fXQp{p) = Xop(p), 

hence we have [l^ 

Xop(p) = (f>Qp{p)uAQ{v,s), 

where 4>Qp{p) is a scalar function. 



(26) 



(27) 



(28) 



However, fxip{p) and fX2p{p) do not equal Xip{p) and X2p{p) respectively, so we define 



xtp{p) = ^-^xipb), 



1 -i> 

Xip{p) = -^xip{p), 

xtpiP) = ^y^X2p(p), 
1-4 

X2P(.P) = —^X2P{P)- 



It is easy to see that 



l^xt2p{p) = Xt2p{p)^ 

i'Xi,2p{p) = -xT,2p{p)- 

Hence, like xop{p), Xi2piP) '^an be expressed as 

Xt2p{p) = '^mp(p)waq(w,s), 

where 4>ip{p) and (j^tpip) are scalar functions. 
By using Eq. (26) and (31), we can see that 

Xl,2pip) = tt4>l,2pi.P)'^f^Q{v, S), 

where again 4>'^p{p) and 4>2p{p) are scalar functions. Then we obtain equations for all the scalar functions: 



^^pW = ^„,^^^^/^„^,vU.2_M/2^,v^ / 7:5^[^3 + (p/ + «)^4]<^op(9) 



{pi +mD + Eo + ie){p'f -W'^+ ie) J {2n} 
—i 

+ 



{pi + mD + Eo+ ie){pf -Wl + ie) 
2(p; + mz)+i;o+«£)2(pf -VK^ + ie) 7 {^-ny 



'-(^^ = 2ip, +mn+Eo +ie){pf - + ie) / (S^^^^ + + m)V2]<t>0P{<l), 
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(j}i + mD + Ea+ie){pf -W^ + ie) J {2tt)^ 
{pi + + i;o + «£)(P; - VVp + K) J (27r)* 



i{-2Er+p^,) 



2{p, + + Eo + ^enpf - W^, + ^e) J (2.)^^''' + + mm^,p{q) 

d*q 

2{pi + niD + Eo + ie){pf - + le) J (2ttP 



-i{-2E^+p1) 



dS 

2{pi+mD + Eo + ief{pf-Wl + ie)J (2^ 
ipl 



^{pi + mo + Eo + ie){p} - Wl + ie) ' {pi + mo + Eq + ie f{pf - + ie) 



[Vs + {pi+qi)V4](l)op{q) 
iE2 



-li2Ei-piy 



d^q 

Mpi +mD + Eo + is)3{pf -Wl + ie)\ J (2^ 
—I f d'^q 



[Vi + [pi + qi)V2](j)ap{q) 



2{pi + niD + So + ie) {pf -Wj, + is) J {2tt) 



j[Pt ■ qt -Pt\Vi(j3ap{q)-, 



(37) 



2(p2 -Wl + ie) 
+ 



'^''^ [Vi + {pi+qi)V2]^C^ip{q) 



{2n 



2{pi + niD + Eo+ ie){pf -W^ + is) 
i f d*q / Vf qt 



Pt -Pt 
d^q 



2{pf-W^ + is)J (27r) 



1 



Pt - Pt 



(2^)4 

V4(l)Qp{q) 
d^q 



[Vi + {pi + qi)V2]^tpil) 



2(P, 



•I + niD + Eo + is)^{pf 
d^q 



- IS] 



(27r) 



[Vz + {pi + qi)Vi]<f)op{q) 



wl 



- IS) J (27r)4 
t{2E,~pj) 



[Vi + [pi + qi)V2\(t>op{q) 
dS 



4{pi +mD + Eo+ is)'^{pf -Wl + is) 
From Eqs. (34) and (36) it is easy to see that [17 1 



(27r) 



j[Vi + {pi+qi)V2\^0P{q). 



^ipb) = \<t>op{p)- 



(38) 



(39) 



The numerical solutions for (j)Qp{p) and (f>^p{p) can be obtained by discretizing the integration region into 
n pieces (with n sufficiently large). In this way, the integral equations become matrix equations and the BS 
scalar functions (popip) and (j)'^p{p) become n dimensional vectors. Thus (f>op{p) is the solution of the eigenvalue 
equation (A — I)(j)Q =0, where A is an n x n matrix corresponding to the right hand side of Eq. (34). In order 
to have a unique solution for the ground state, the rank of {A — I) should be n — 1. From Eq. (35), (f)'^p{p) is 
the solution of (A — I)4>i = B, where B is an n dimensional vector corresponding to the first and third integral 
terms on the right hand side of Eq. (35). In order to have solutions for (p^p{p), the rank of the augmented 
matrix {A — /, B) should be equal to that of {A~ I), i.e., B can be expressed as linear combination of the n — 1 



9 



linearly independent columns in {A — I). This is difficult to guarantee ii B ^ 0, since the way to divide {A — I) 
into n columns is arbitrary. Therefore, following Ref. pLZj we demand the following condition in order to have 
solutions for (j)lp{p) 



Vs. + {pi + qi)Vi + 



pi + ruD + Eq - 
The simplest forms for V3 and V4 satisfiying Eq. (40) are 



■ le 



pi + WLD + Eo + ie 



Vi + ipi+qi)V2 



(40) 



(41) 



V4 - 



(42) 



pi + mo + i?o + 

With Eq. (40), (j>^p{p) satisfies the same eigenvalue equation as (t)Qp{p). Therefore, we have 

(t)+p{p) = cT<Pop{p). (43) 
where cr is a constant of proportionality with mass dimension. It was determined to be zero by applying Luke's 



theorem 



15| at the zero-recoil point in HQET 17|. Therefore, (f)^p{p) does not contribute 



As to the forms of and Vg , we simply assume that 

^6 = 5''V2. 



(44) 



where (5 is a parameter which is of order Aqcd- 

From Eq. (38) we can see that (p^pip) is only related to 4>op{p) and 4''^p{p), so our goal is to solve the BS 
equations of (f>op{p) and (f)2p{p) numerically. 

In the covariant instantaneous approximation, = ^^1^,=^, , i = 1, 2 [l^-[3|. In the meson case, one has [3| 

d^k Sttk' 



Sttk' 
[{Pt - qtY + 



212 



(27r)3 (fc2 + //2)2' 



167r 



a 



2 \ meson 



(45) 



(46) 



3 {pt - qtY + fJ.'^ ' 

where k' and ctseff are coupling parameters related to the scalar confinement and the one-gluon-exchange terms, 
respectively. The parameter is introduced to avoid the infrared divergence in numerical calculations, and in 
the end we will take the limit /i — > 0. In the baryon case, since scalar confinement term is still due to scalar 
interaction, the form of Vi need not be changed. Only the parameter k' in the meson case has to be replaced 
by K which describes the scalar confinement interaction between the heavy quark and the diquark. However, 
the diquark is not a pointlike object, there should be a form factor in V2, F{Q^){Q = p2 —^2)1 to describe the 
structure of the diquark (see FIG.l) [l9|. 



2 ' 



(47) 



where Qq is a parameter which freezes F{Q'^) when is very small. In the high energy region the form 



factor is proportional to which is consistent with perturbative QCD calculations 20 1. By analyzing the 
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TABLE I: Values of ctsej j for different mo and hi. 

K{GeV^) 0.02 0.04 0.06 0.08 
mi3 = 650MeF a^e// 0.66 0.7 0.73 0.75 
mo = 700MeV a^e// 0.69 0.73 0.76 0.78 
mo = 750MeV a^e// 0.74 0.77 0.79 0.8 
mi3 = 800MeV' a^^// 0.77 0.8 0.82 0.83 



electromagnetic form factor for the proton, it was found that Qq = 3.2GeF^ can lead to consistent results with 
the experimental data [l9j. Based on the above analysis , the form of kernel for BS equation in the baryon case 
is taken as 161 



Sttk 



212 



(27r)3 (fc2 +Ai2)2^ 



(48) 



V, = -1^ "-//Q" (49) 

From the BS equation solutions in the meson case it was found that the values mi, = 5.02GeV and to^ = 
1.58GeV give predictions which are in good agreement with experiments [iS]. Hence in the baryon case we take 

mD+Eo + —El + -^E2 = 0.62GeK (50) 

nib ml 

The dimension of k, is three and that of k' is two. This extra dimension in k should be caused by nonpertur- 
bative diagrams which include the frozen form factor F{Q^) at low momentum region. Since Aqcd is the only 
parameter which is related to confinement, we expect that 

K = Sk', (51) 

where S is of order Aqcd- It is noted that the two proportionality parameters in Eqs. (44) and (51) could be 
different in general although they are both of order Aqcd- However, in order to simplify our model we assume 
these two parameters are the same. The value of k was determined to be in the region between 0.02GeV^ and 
O.OSGeV^ in Ref. [IJ. Since k' is about 0.2GeV'^ we take S to be in the range from O.lGeV to OAGeV. 

One notes that Ei ^ AqcdEq, E2 ^ AqcdEi, hence we further assume that Ei — 6Eq, E2 = 5Ei = S^Eq. 

In general, ^op(p) can be the function of pi and pt- Defining (j)Qp{pt) = J{dpi/2T:)(f>Qp, one gets immediately 
the BS equation for (pQp{pt)- 



1 f d^qt 

2{-Wp + fUD + Eo)Wp J (27r)3 



^opiPt) = - o^ w , , z.„M.. / T^Jb (^1 ~ 2WpV2)Mqt), (52) 



where we have used the residue theorem and selected the upper contour in the pi plane which has a singular 
point -Wp + ie (see FIG.2). 
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FIG. 2: Three singular points, —mo — Eq — ie, —Wp + ie and Wp — ie, in the pi plane. 



Substituting Vi and V2 into Eq. (52) we have 

{-Wp + niD + E^)(j)Qp{pt) = - 



1 I f qfdqt 



2W^ 



^op(gt) + 



3(Oo-u2) 



gf^gggf 1 
47r2 2ptqt 

1 /■ 9?c^<7t 



'{pt - gt)2 + u2 



47r2 (p2 +q2 +u2)2_ 4p2^2 



4>op{qt 



(53) 



This is a eigenvalue equation of (j)op{pt)- Solving this equation we obtain the values of the model parameter 
aseff corresponding to different values of k and rriD (see TABLE 1). 

With the numerical result of 4'op{pt), we can solve out (f)2p{pt) numerically from Eq. (37). 
The normalization equation for the BS equation takes the following form 21|: 



(27r) 



d\dWxk{q')4riI{q\ q, k) + K{q' , q, k)]xk{q) = 2fco, 
Okn 



(54) 



where I{q' ,q, P) = 6^^''{q' — q)[SFipi)]^^[SDi—P2)]^^, in which pi = \iP + q' ,p2 = —X2P + q', and feg represents 
the energy of Ab, which is the mass of Af, in the rest frame of A;,. With this condition we can give the plots of 
normalized <j)op{pt) and (f>2p{pt) in FIGs. 3 and 4, respectively. 




P, (GeV) 



FIG. 3: Plot of 4>opipt) with mo = 650MeV, k = 0.04. 
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FIG. 4: Plot of (?!)+, (pt) with mo ^ 650MeV , k = 0.04. 



VI. Numerical results of ej, 

With the numerical results of the BS amplitudes solved out above, we can calculate the value of eb- The left 
hand side of Eq. (5) can be expressed as the overlap integral of the BS amplitude as the following: 



Hence we have 



{Ab{v,s)\b-f^j5b\Mv,s)) 



(1 + ef,)uA, (v, s)7-^75UA(, (v, s) 



(27r) 



^Xp{p)i^15Xp{p)Sd{P2)- 



■Xp{p)l^lbXp{p)S jj^ {P2)- 



(27:)' 



Expanding both sides of Eq. (56) to order , we have 



rrib 



_dS>_ 
(2^ 



^^P^P) + -^i^lpiP) + <t>ip{p)tt) + ~2 (<^2P(P) + <l^2p{p)Pt 



X7N5 



4>op{p) + —{4>^p{p) + <P^p{p)f^) + \Utp{p)<p2pip)ft 



UA^iv,s)Sjj^{p2), 



where we have expanded eb in Eq. (5) as 



(55) 



(56) 



(57) 



Ef) = £06 H Clb H o 

rrib mf 



£26- 



(58) 



By comparing the two sides of Eq. (57) at each order in the l/niq expansion, we have the following equations: 
to leading order in the 1 / mq expansion: 



(1 + eob)u\^ (w, s)7^75'"Ab (w, s) 



d^p 



UA^V, s)(pop{p)'y^j5(l)op{p)uA,{v, s)Sp,^ {P2) 



(59) 
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TABLE II: Values of £26 for different mo and k. 

ii{GeV^) 0.02 0.04 0.06 0.08 

mo = 650MeV e2b{GeV^) -0.11 -0.19 -0.26 -0.32 

niD = 700MeV e2b{GeV^) -0.16 -0.25 -0.32 -0.41 

niD = 750MeV e2b{GeV'^) -0.27 -0.36 -0.44 -0.49 

ruD = SOOAfeV e2b{GeV^) -0.39 -0.47 -0.56 -0.62 

to first order in the I/toq expansion: 



+</'ip(p)7^75'/'op(p) + 0rp(pVt7^750op(p)]wAi,(w,s)S'^^(P2), (60) 



to second order in the l/mg expansion: 



e2bMA,(w, s)7 75MA,(i;,s) = y ^^;^WA,(w, s)[(/>op(p)7 750Ip(p) + <?!'op(p)7 75^j02p(p) + </'rp(p)7 750rp(p) 

+<^]*p(p)7'^75|>t0rp(p) + 0rp(p)|*t7'^750jp(p) + 0];p(pVt7'^75|>j0rp(p) 

+'^2pb)7'^75'/>op(p) + 02pb)|>f7^75'/'op(p)]wAi,(w, s)5^^(p2)- (61) 

Getting rid of the spinors, ef, at each order in the l/rriQ expansion can be expressed as the foUowing: 

(1 + eob) - J ^Mp)Mp)SD'iP2), (62) 

j^^(t>op{p)(t)tpip)SD^{p2), (63) 

= 2 [ ^Mp)^tp{p)(.Pi ~Wl)~\( ^J,-p[p)^-p{p)pl{pf - Wl). (64) 



From Eq. (62) one can see immediately that eob = because the Isgur-Wise function is normalized to 1. From 



Eq. (63) we have ei^ — since (l)fp{p) = 0. This is consistent with the demand of the current conservation [15j. 

Substituting Eqs. (34), (37), and (39) into Eq. (64), using 0(pt) = / ^4){p) and the residue theorem, we 
obtain the value of e2b with the obtained numerical results of <j>ap{pt) and (p2p{pt)- The results are listed in 
TABLE 11. 

We can see from TABLE II that the value of £26 depends on our model parameters, mu and k. In the ranges 
of the parameters, e2b is always negative. The absolute value of £26 increases with the increase of rriu and k. 
When mD = 650MeV, ^Ie^ changes from OMGeV to 0.56GeV. When mo = SOQMeV, y^\e^\ changes from 
0.62GeV to 0.796*6^. As expected, ^/\e^\ is of order Aqcpi- This indicates that our results are reasonable. 

IV. Summary and discussion 

More and more data have been and will be collected at LHC. This makes more precise measurement of 
Cabibbo-Kobayashi-Maskawa matrix elements, especially the value of Vub, become possible. Vub can be measured 
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through the inclusive semileptonic decays of polarized Af,. Since the decay rate involves two parameters, /Li^ 
and Cb, to second order in the l/rub expansion, the theoretical calculation of these two parameters is important 
to extract more precise Vub from experimental data. /LtJ was calculated in the previous work. In this paper we 
focus on the theoretical calculation of £(, in the BS equation approach. Since £(, is only non-zero at second order 
in the 1/mb expansion, it is necessary to establish the BS equation for to this order. 

We regard the heavy baryon A^ as composed of a heavy b quark and a light scalar diquark based on the 
fact that the light degrees of freedom in A;, have fixed spin and isospin quantum numbers. In this picture, 
we establish the BS equation for A^ to second order in the 1/mb expansion. Then we solve the BS equation 
numerically by applying the kernel which include the scalar confinement and the one-gluon-exchange terms. To 
second order in the l/rrib expansion, we obtain numerical results for BS scalar functions. Expressing as the 
overlap integral of the BS amplitude for A^ we obtain the numerical values for £{,. It is found that is only 
non-zero at l/m^ as expected. This indicates that our result is consistent with the requirements from the heavy 
quark symmetry and the current conservation. At eb can be expressed as e2b/m^. We found that in the 

ranges of our model parameters, €2b is always negative, and varies form —O.llGeV^ to —Q.62GeV^ . We find 



that \/\e2b\ is of order Aqc d as expected. 

There are some uncertainties in our model. Compared with the heavy meson case, heavy baryons are much 
more complicated. To order 1/mb there are two parameters in our model, m^ and k 17]. When we expand to 
second order in the 1 / mt, expansion, more uncertainties are involved, including the forms for the kernel V5 and 
Ve and the values of binding energies at orders 1/mb (Ei) and (£^2)- To simplify our model we assume 

that V5 = (5^Vi, Vg — 6^V2, E2 = 5Ei = S'^Eq, with the same proportionality parameter 6 which is of order 
^QCD- Our model to first order in the l/mq expansion has been tested in A;, — > Ac semileptonic decay and 
Ab — ^ AcTT nonleptonic decay 0[l2| 22 1. Further tests are needed in other physical processes and for the model 
to l/m2 . 
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